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In this paper, a generalization of standard spin fluctuation theory is considered 
by replacing the simple Hubbard interaction by the screened Hartree-Fock inter- 
action for f electrons. This model is then used in both an LS and a JJ coupling 
scheme to construct the particle-particle scattering vertex in an on-site approxima- 
tion. This vertex is shown to lead to an instability for a superconducting pair state 
which obeys Hunds rules, with L=5, S=l, and J=4. The degeneracy of this state 
is broken by anisotropy of the quasiparticle wavefunctions. Detailed calculations 
are presented for the case of UPt%. 
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After over a decade's worth of theoretical work, there is no overall agreement on a 
microscopic theory for heavy fermion superconductivity. The overall prejudice, though, 
is that the underlying pairing mechanism is similar to that operative in superfluid ^He. 
Most attempts at a theory based on this approach have been to make the simplest possible 
modifications to the standard single orbital Hubbard interaction used in the He problem. 
These attempts have had mixed success. The philosophy of this paper will be to actually 
do for heavy fermions what was done for 3 ife, that is to use the Hartree-Fock interaction 
between f electrons including full orbital and spin-orbital effects and construct the effective 
particle-particle vertex by the appropriate diagrammatic summation. In principle, this 
theory contains all relevant physics within a spin fluctuation based approach. Even at 
the simplest level, new physics emerges which is not present when using a simple Hubbard 
interaction. In particular, the maximum instability in the particle-particle vertex occurs for 
a pair state which obeys Hunds rules. For f electrons, this corresponds to a state which has 
L=5, S=l, and J— 4. The degeneracy of this multiplet is broken in real metals by crystalline 
anisotropy effects in the normal state. This is reflected by (1) the orbital and momentum 
dependence of the bare susceptibility bubble which forms the internal lines of the vertex 
and (2) the orbital and momentum dependence of the quasiparticle wavefunctions which 
form the external lines. In this paper, (1) is treated in a simple manner and (2) is treated 
within a band theoretic approximation. The frequency dependence of (1) and (2), which 
acts to set the overall scale for T c , is also treated in a simple fashion. The theory has the 
advantage that it can be systematically improved by removing these approximations. The 
above ideas are illustrated by calculations for UPt^. 

In the first section, a motivation of this theory is given by looking at some systematics 
of heavy fermion superconductors and by comparing the heavy fermion problem to that 
of He. In the second section, the general formalism is described. The single orbital 
version of this theory is shown to yield the paramagnon model for 3 iiZe. The standard spin 



3 

fluctuation models worked on previously are then shown to be a lattice generalization of 
the single orbital model. In the third section, the formalism for the f electron problem is 
derived, with the particle-particle vertex equations solved for various approximations for 
the susceptibility bubble in both an LS and a JJ coupling scheme. General properties of 
the vertex are described based on group theory. In the fourth section, the pair vertex is 
projected onto the Fermi surface. Calculations are then described for the case of UPt% 
utilizing information from a relativistic band structure calculation. In the last section, 
future directions, including the question of inter-site pairing effects, will be discussed. A 
shorter version of this work has appeared earlier. 1 

I. Introduction 

Sufficient evidence has accumulated over the past eleven years to demonstrate that 
the superconductivity seen in a number of f electron metals with large effective mass is 
unconventional in nature, that is, the group representation describing the order parame- 
ter is almost certainly not the identity representation (i^). This, along with a variety 
of other facts, casts doubt on a traditional electron-phonon mechanism as mediating the 
pairing. The first theoretical work in this area ten years ago showed a close connection 
of these metals with superfluid 3 .He. 2 In particular, they are near both a magnetic and a 
localization instability. A classic example is UPt%. Doping with Pd, for instance, causes 
this metal to become strongly antiferromagnetic. Further doping causes the f electrons 
to become localized. 3 Anderson 2 also emphasized that the on-site part of the interaction 
must be playing a major role given the large ratio (~0.1) of the superconducting transition 
temperature to the Fermi energy. This important observation has been largely forgotten. 
Anderson 4 was also the first to point out that heavy fermion superconductors have two f 
atoms per unit cell. For on-site pairing, one can have an odd parity ground state in this 
case (with one atom per cell, the pair state would have to be even for on-site pairing). 
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This unusual observation has also been largely forgotten, except in a later paper by Appel 
and Hertel 5 where a formalism for describing localized pairs for UPt^ was developed in 
great detail. The reason the above points were largely forgotten was the observation of 
antiferromagnetic spin fluctuations in several heavy fermion superconductors by neutron 
scattering.^ Such fluctuations occur because of exchange interactions between near neigh- 
bor sites. This led to a picture of near neighbor pairing based on these fluctuations by a 
number of authors. 7 ' 8 Subsequent work largely concentrated on generalizing these simple 
theories to handle the non-symmorphic (HCP) lattice structure of UPt%?^ These theo- 
ries have had mixed success. In particular, available data on U Pt% point to a pair state 
from a two dimensional group representation with both line and point nodes and proba- 
bly of odd parity. ' This state would have r§ (i^u) symmetry. Although non-trivial 
group representations occur in these calculations, this particular state has never emerged 
as the ground state. The last work done in this area by the author 1 ^ indicated that the 
anisotropy of the quasiparticle wavefunctions plays a fundamental role because this prob- 
lem cannot be reduced to an effective one-band form given the two f atoms per unit cell. 
Therefore, simple theories as pursued above will simply be inadequate for describing real 
heavy fermion metals and any results generated by them questionable. 

This paper is an attempt to break this theoretical deadlock. This work was motivated 
by several additional issues than those listed above. A number of alternate theories have 
been proposed recently, in particular by Cox 14 and by Coleman et al, 15 which emphasize 
an on-site pairing viewpoint. Cox's work is important in that he emphasized the important 
role that orbital effects play in this problem. Another key motivation was an experiment 
by Osborn et al 16 which detected excitations between Coulomb multiplets with high en- 
ergy neutron scattering, not only in localized f metals like Pr and UPd%, but also in 
U Pt% itself. This indicates that multiplet correlations present in atoms survives even in a 
metal with itinerant f quasiparticles. In Table 1, a list of the seven known heavy fermion 
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superconductors are listed. There are two striking things about this table. First, six of 
the seven known heavy fermion superconductors are uranium alloys. Moreover, there is 
strong experimental evidence that the uranium atoms are close to an f 2 configuration. 
The magnetic susceptibilities of UPt^ and UPd^Al^ 1 look almost identical to that of 
the local / metal PrNi§. The susceptibility of URU2S12 has been most successfully 
explained based on an f 2 ground state. 19 Cox's quadrupolar model for UBe\^ 2Q is also 
based on an f 2 configuration. It should also be noted that UPt^ is very similar to UPd^ 
(similar crystal structures, almost identical f atom separations) yet the latter is clearly a 
local f 2 metal. 21 In the high energy neutron data, 16 the Coulomb excitation seen in these 
two metals looks very similar. This would be hard to accept if UPt^ was not close to being 
f 2 . As for CeCu 2 Si 2 , it may not be like the rest, although it has been pointed out that 
an f 2 admixture is needed to explain its properties with the Anderson impurity model. 22 
The importance of these facts is that since the f atom has two bare f electrons per site, 
this leads to a strong motivation that the superconducting pairs have two f quasiparticles 
per site from a trial wavefunction point of view. The second striking point of Table 1 is 
that all of these metals either have two f atoms per unit cell, or undergo some sort of 
magnetic transition at temperatures above the superconducting transition which gives a 
new unit cell with two f atoms. As discussed above, this fact has little relevance to a near 
neighbor pairing model (since an atom in any crystal structure has near neighbors), but 
plays a crucial role for on-site pairing (since one can have even parity or odd parity pairing 
depending on the phase of the order parameter on the two sites) . 

A further motivation of the importance of on-site pairing can be obtained by comparing 
the case of 3 He to uranium alloys. In Figure 1, a plot is shown of the interaction potential 
of two He atoms. 23 This is amazingly similar to what one would expect of two f electrons on 
a uranium site. In particular, there is strong repulsion at small interparticle separation due 
to the Coulomb repulsion between the two f electrons, there is attraction at intermediate 
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distances (or order 3 a.u.) since the ion core attraction exceeds this repulsion in this range 
(which leads to an / 2 ground state), and then the potential weakens at large separation due 
to the exponential decay of the f electron radial function. The direct interaction potential 
for 3 ife, though, has been shown to be inadequate for describing the pairing of He atoms 
in the superfluid state (it predicts L=2 pairing). 24 The reason for this is the major role 
that collective effects play due to polarization of the medium. 25 ' 26 ' 27 ' 28 ' 29 This led to 
the development of a paramagnon model for 3 ife. In this model, a much simpler direct 
interaction potential is used, a repulsive contact interaction between atoms of opposite 
spin. But this potential in turn is used to sum a diagrammatic series to all orders, thus 
including the important collective effects. This gives a good description of the superfluid 
state of 3 i7e. 27 ' 31 - 1 In the next section, a generalization of this model is considered which 
includes orbital interaction effects necessary in dealing with f electrons. It is then shown 
that the paramagnon model for 3 ife, as well as previous spin fluctuation models for heavy 
fermions, are subsets of this more general theory. 

II. General Formalism 

The particle-particle vertex is defined by 

pabcd _ pabcd _ paecf _^ef pfbed ^ 
e,f 

where i~b is the bare vertex, xo is the bare susceptibility bubble, the indices label orbitals, 
and the minus sign is due to the closed Fermion loop defining the bubble (in T, the 
first two indices label incoming lines, the last two outgoing lines). i~o is taken to be the 
antisymmetrized Coulomb interaction (y a ^ cd — V abdc ) 

r^abcd \ abed rp /o\ 
7 - 2_^ c k t k i 2 ) 

k 

where cy. are combinations of 3j symbols and are Coulomb multipole (Slater) integrals 
defined on page 217 of Condon and Odabasi. 31 (A simpler expression of this type has been 
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used in earlier spin fluctuation work. 32 ) For s electrons, Eq. 2 reduces to 

r bcd = Sachd ~ S adhc = ^ac?>bd -°ac-o U ) (3) 

where the last term is a scalar product of Pauli spin matrices and the indices now label 
just spins (f for up, 2 for down). Eq. 1 is easily solved, giving 27 ' 28 

r 1111 = -f 2 xo/(i - F$xl) (4) 

r 1212 = Fo + F 3 X 2 /(l-F 2 ^) (5) 

Eq. 4 (5) is a sum of odd (even) number of longitudinal (e=f) bubbles. The triplet (S=l) 
vertex is just Eq. 4, the singlet (S=0) vertex is 2P 1212 - T 1111 (since T 1212 is half the 
sum of the singlet and triplet vertex 27 ). Alternately, the singlet vertex is pl^ 2 — pl221 
(i.e., antisymmetrizing Eq. 5) where the latter term is a sum of transverse (e^f) bubbles 
(i.e., ladder diagrams). These expressions are equivalent of course. Eqs. 4 and 5 form the 
basis of the standard paramagnon model (where Fq is generally denoted as I). 

The first thing to note from the above is that the triplet vertex is negative (attractive) 
and the singlet one positive (repulsive) . Thus, the vertex exhibits Hunds first rule (maximal 
S). In the 3 iiZe problem, the L of the pair state is determined by projecting Eq. 4 onto 
the Fermi surface. Since it is isotropic, and L must be odd (since S=l), it is necessary to 
include the momentum dependence of xo to obtain a non-zero projection. As expected, 
L=l is found since this function has the largest projection on a spherical Fermi surface for 
odd L harmonics. Note that the momentum dependence is not critical, it is only neccessary 
to give a non-zero projection. In fact, very different models of the momentum dependence 
of the vertex give identical pairing coupling constants. 33 A physical picture of the pairing 
in these models based on mutual interaction of the two particles via their polarization 
clouds has been given by Leggett. 29 

This picture can be contrasted with that given for the heavy fermions and high T c 
cuprates based on singlet (L=2) pairing. 7 This "violation" of Hunds rules is obtained 
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by considering pairing of electrons on near neighbor atoms in a nearly antiferromagnetic 
metal. In such a case, x(Q) > x(0) where Q is the ordering wavevector (assumed to 
be commensurate with the lattice corresponding in real space to anti-alignment of near 
neighbor spins) and x is the dressed bubble. This means in real space that x(-R, R') is 
negative (R,R' are site indices, with R a near neighbor of R). This can be achieved by 
having momentum dependence in either xo or in i~b (x being defined by an equation similar 
to Eq. 1). The latter is preferable, in that commensurate Q are rarely obtained for xo 
except in special circumstances, and has been used for fitting neutron scattering data both 
in heavy fermions 6 and in high T c cuprates. 34 Now, the same thing which gives a negative 
x{R, R') also gives a sign for r{R, R' \ R, R') opposite to that of -T(-R, -R, R, R). Therefore, 
if the order parameter is such that A(R, R) is zero and A(R, R') is non-zero (such as found 
for certain d-wave states), then one can have S=0 pairing. Note the complete difference 
in physics than that discussed above for the ^He problem. In particular, the momentum 
dependence is crucial for this argument. From the lattice point of view, ^He is actually 
more closely related to on-site pairing models than to near neighbor pairing models. 

III. f Electrons 

The formalism of the previous section is now applied to the problem of f electrons. 
At the bare interaction level, we are already faced with the problem that the uranium 
ion is in the intermediate coupling regime (i.e., midway between LS and J J coupling). 35 
On the other hand, the spin-orbit interaction is large enough that in electronic structure 
calculations, no j=7/2 quasiparticles are occupied. 36 Because of this, the susceptibility 
bubble xo wm be almost pure j=5/2 in character. Therefore, even if the bare interaction 
is in the intermediate coupling regime, the effective interaction for quasiparticles which 
comes out of Eq. 1 will be in the JJ coupling limit. Despite this, we will start out by 
deriving results in the LS coupling limit to make connections to the 3 ife problem discussed 
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in the last section. Then, we will turn to the JJ scheme. 

We start by reviewing the multiplet structure of the / uranium ion, shown in Figure 
2.3^ There are 3 spin triplets (L=l,3,5, each spin-orbit split into 3 orbital multiplets) and 
4 spin singlets (L=0,2,4,6). This level structure can be fit by the following scheme. At 
the Hartree-Fock level, only Slater integrals of even rank appear. Fits of the spectra can 
be achieved by reducing these integrals on average by 38%. This effect is due to screening 
caused by Coulomb correlations (i.e., the particles try to avoid one another, thus reducing 
their effective interaction) which can be approximately calculated within a configurational 
interaction (CI) scheme. In addition, CI causes effective operators of odd rank to appear 
not present at the Hartree-Fock level (known as Trees parameters). These terms are rather 
small, though, and we ignore them. To discuss the level scheme, it is useful to find linear 
combinations of the coefficients of Eq. 2 which more clearly reflect the group theoretical 
structure of the f electrons, which are labeled e^. This has been achieved by Racah 37,38 ' 31 
and is equivalent to replacing the Slater integrals Fj~ (k=0,2,4,6) by linear combinations 
Efc (k=0,l,2,3). Eq is defined such that all f 2 terms have this energy (eo=l). It is equal 
to Fq plus a linear combination of the other terms and is equivalent to the Hubbard U 
parameter (the other E^ parameters do not contain Fq). E\ is defined so as to distinguish 
spin singlets from spin triplets, with ei=0 for triplets and ei=2 for all singlets but L=0. 
Note that these coefficients are identical to the s-electron problem of the previous section. 
Thus, E\ plays the same role in the f electron problem as the paramagnon I (Fq) plays in 
the 3 _£fe problem. This is quite interesting, since E\ is a shape fluctuation term instead 
of a charge fluctuation term (i.e., it does not involve Fq). This means that the inclusion 
of orbital degeneracy qualitatively changes the physics relative to single orbital models. 
As for the L=0 singlet, it has e\=9. This is another consequence of orbital degeneracy, 
basically L=0 has a different group structure than L=2,4,6 since it already appears at the 
/° level (i.e., it has a different "quasi-spin" or "seniority"). We summarize by writing 
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down the expression for e\ 

ei =< LS\q 12 + X - - 2si • s 2 |LS >= 2 - 5(5 + 1) + 7S L0 (6) 

where 1,2 label the two electrons, q is the seniority operator, and s is the spin operator. 
Note the similarity to the last expression of Eq. 3. As for e 2 , it is isomorphic to L=2,4,6 
and thus acts to split these three multiplets apart (it is zero for all other L states). The 
expression for this term is quite complicated and will not be written down. The most 
interesting parameter is E%. It acts to split apart the three spin triplets, with e3=-9,0,33 
for L=5,3,l respectively (it is non-zero for all terms but L=0). This can be written as 

e 3 = [±-S(S+l)][L(L+l)-24g(U)} (7) 

where g(U) is the Casimir operator of the group G 2 with U labeling the representation 
of G 2 appropriate for a particular LS state (note that L(L+1) is the Casimir operator for 
the group SO3). The interest is that this has similarities to the case of p electrons, where 
VanVleck 39 showed long ago that the interaction between two p electrons can be written 

as 

W12 = F + (-5 - 3h ■ h - 12si • s 2 )F 2 (8) 

— * 

with / the orbital angular momentum operator. This looks like an orbital generalization 
of the last expression in Eq. 3. A summary of the coefficients for the / 2 states is given 
in Table 2. 

We now turn to a solution of Eq. 1. At this time, we will assume that xj = xo (i- e - 5 n ° 
orbital dependence to the bubble). This approximation will be discussed below. States of 
definite LS have antisymmetrized wavefunctions. So, the vertex for the state L=5, Mi=b 
(S=l, M s =l) will be 

7-1 7-i3232 7-i3223 7-i3232 /r\\ 

where indices label m; (with all spins up). The first term involves longitudinal bubbles, 
the second transverse bubbles, and the last compensates for double counting (since the 
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bare vertex is antisymmetric by definition). Vertices for other LS states can be obtained 
by either using the appropriate antisymmetric combination of the r or by employing 
Slater's diagonal sum rule. 31 This is equivalent to the two forms of the singlet vertex in 
the paramagnon model discussed after Eq. 5. 

In certain cases, analytic results can be derived by expanding Eq. 1 into a bubble 
summation, just as was done in the previous section for 3 ife. In particular, let us start 
with just including the Eq term. The result for all states of definite LS is 

r = E /[(i - e 0X q)(i + uEqxo)] + £ 2 xo/(i - ^oxo) (10) 

The first term comes from longitudinal bubbles, the second from transverse bubbles and 
double counting. Note that 13 is the orbital degeneracy (14) minus 1. This expression, 
plotted in Figure 3, is always repulsive. The behavior of this term is that as xo increases 
from zero, the repulsion is reduced compared to Eq then begins to increase again and 
diverges at Eqxq=1- This corresponds to a localization instability. 

Now assume that only E\ is non-zero. Analytic results can also be obtained. For the 
spin triplet states, Eq. 1 is now 

r = -ll£?fxo/[(l " 81Efxo)(l - ±Elxl)} + 2Efxo/(l ~ ^ X q) (11) 

where the first term comes from longitudinal bubbles, the second from transverse bubbles. 
This expression, plotted in Figure 4, is zero for xo=0 an d then has a negative divergence as 
9E\xq approaches 1, corresponding to a magnetic instability. This behavior is analogous 
to the triplet vertex in 3 i/e. Eq. 1 can also be solved for the L=6 singlet 

r = {AE\ + 13^X0 " 126£fxo " 162£i 4 Xo)/[(l - 81£?Xo)(l - ^xq)] ~ 2E 1 (12) 

where the first term comes from longitudinal bubbles, the second from transverse bubbles, 
and the third from double counting. Just as for the singlet vertex in 3 ife, this interaction, 
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plotted in Figure 4, is repulsive for xo=0 (equal to 2E\) and has a positive divergence at 
9Eixo=l- This expression should also be valid for the L=2 and L=4 singlets. The vertex 
for the L=0 state would be more repulsive. 

The general solution to Eq. 1 is difficult to construct analytically due to the compli- 
cated nature of the bare vertex when all terms are included. Solving numerically would 
also appear to be difficult since their are four orbital indices involved. Progress can be 
made, though, by defining 

pbe pb+n,e,b,e+n (13) 

with this definition, Eq. 1 reduces to 

r be = rS e -J2r b f xl J+n r^ (u) 
f 

Thus, an iV 4 matrix equation has been reduced to an iV 2 matrix equation for each value of 
n (n=0 are the longitudinal bubbles, the rest are transverse bubbles). This is easily solved 
on the computer given input values for the E^ and xo- 

The Ejg parameters were taken from Goldschmidt. 35 They have been fit to uranium ion 
data (the level scheme for xo=0 in Figure 5 differs in some quantitative details from the 
experimental level scheme of Figure 2 since spin-orbit has not been included at this point). 
These parameters are on average 62% of their Hartree-Fock values due to screening. By 
comparing to the high energy neutron scattering data of Osborn et al on UPt^, 1 ® these 
parameters should be reduced by another 28% when going into the solid. The latter effect 
is only of quantitative significance and is based on one transition which is seen (assumed 
to be from 3 i?4 to 3 -F2), so is ignored in order to use the well established values of E], 
listed in Goldschmidt's article. Hopefully, detailed solid state values of these parameters 
will become available with further experimental work. An additional note is that the fitted 
Eq is referenced to some arbitrary value of the energy zero, so has no intrinsic meaning. 
The value listed in the Goldschmidt article, though, gives an Fq of 1.83 eV which is 
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fortuitously close to estimates of the screened Coulomb U for uranium, 40 so we retain it 
without adjustment. 

A final comment concerns the energy zero of the problem. Superconductivity involves 
an instability of the Fermi surface. For a uranium ion, two f electrons are occupied. Even 
in band structure calculations for UPt^, the number of occupied j=5/2 f electrons is just 
above two. Thus, the term Eq (the Coulomb repulsion between the two f electrons) is 
already included in the definition of Fermi energy and represents the zero of energy for 
the uranium problem (for the cerium case, where only one f electron is occupied, Eq is 
not included in the energy zero, since it represents the energy of f 2 above f 1 ). But, the 
energy term for the L=5 ground state of an / 2 ion is Eq — 9E%. What about the term E%! 
Since this term cannot be reduced to an effective single particle form, it would not seem 
to enter into the definition of the quasiparticle Fermi energy. This gives the rather bizarre 
result that the L=5 vertex is already attractive at the bare interaction level. Of course, 
one could imagine a scenario where one considered an / 2 -/ 3 Anderson lattice model, with 
the effects of E% built into the ground states. The effective quasiparticle operators in this 
case might implicitly contain the effects of E$. Since a detailed theory of this has not been 
worked out yet, we cannot make any definitive conclusions one way or the other. Since 
the effective interaction, though, strongly departs from the bare value as xo increases from 
zero, this question is of minor significance. For purposes of this paper, we assume that the 
energy zero is at Eq for the / 2 case and for the J 1 case. 

In Figure 5, r is plotted for the triplet states L=5, L=3, and L=l and for the singlet 
L=6. The other singlets will have similar behavior. Just as found for the 3 .He problem, the 
triplet interactions are attractive and the singlet ones repulsive, with an instability in both 
cases at (Eq + 9.Ei)xo = l- Note this criteria is a combination of the two analytic results 
discussed above, thus the instability has both a localization and a magnetic component. 
This observation indicates that the debate concerning both 3 He and heavy fermions about 
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whether the physics is nearly localized or nearly magnetic is merely semantics, as both 
effects are intertwined. A significant difference from the 3 ife case is the effect of orbital 
interactions in the current problem. In 3 ife, the orbital degeneracy of the pair state is 
lifted by Fermi surface projection of the vertex; in the f electron case, this degeneracy is 
already lifted by the interaction itself. Note that the largest attractive instability is for 
L=5, S=l (in such a state, the Coulomb repulsion is minimized). Thus, the pair state is 
predicted to satisfy both Hunds first and second rules, and is a generalization of the results 
obtained for 3 i/e. 

We now turn to a discussion of the problem at the JJ coupling level, which as argued 
above, is more physically relevant at the quasiparticle level than the LS scheme (or even 
the intermediate coupling scheme). The bare interaction vertex can be gotten by replacing 
the orbital-spin indices of Eq. 2 by the indices \i which range from -5/2 to 5/2 (we assume 
only j=5/2 quasiparticles are involved). By taking into account the mixed spinor nature of 
the relativistic orbitals, the formulas for can be calculated in a manner similar to that 
used on p. 217 of Condon and Odabasi. 31 Note that only k=2 and 4 are involved (k=6 
comes in when considering j=7/2 states). The resulting expressions were checked against 
tabulated results on page 560 of de-Shalit and Talmi. 41 As in the LS case, the are not 
useful in exploiting the group properties of the f electrons, so one rotates to another basis 
Efc. It should be noted that these are not the same Ej, as in the LS case. They have 
been discussed in the context of nuclear physics, where JJ coupling has been traditionally 
of more use, 42 and the analogous coefficients are listed in Table 3. In this case, there 
are only 3 terms, J=4,2,0 (corresponding to L=5,3,l in the LS case). As in the LS case, eg 
is 1 for all states. e\ is 3 for J=0, otherwise which has analogies to e\ in the LS case. In 
particular, J=0 appears at the /° level and thus has a different quasi-spin (i.e., seniority) 
than J=2,4. The term splits the latter two states apart (similar to e$ in the LS case). 
In this notation, the ground state J=4 term is Eq — 5£ , 2 • This is lower than the Eq energy 
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zero just as found for the L=5 LS case. By fitting these values for states of definite J, we 
can infer an interaction vertex analogous to VanVleck's 39 of Eq. 8 for the p electron LS 

case 

wi2 = E + ^E 1 + [-2ji • j 2 - |(1 + 912)1^2 (15) 

where q\2 is a seniority operator (< J\q\2\J >= 65 jq) and j is the total angular momentum 
operator. 

Analogous analytic series can also be constructed. In particular, keeping just Eq gives 
an expression like Eq. 10 with 13 replaced by 5 since the orbital degeneracy is now 6 
instead of 14. 

We now solve Eq. 1 exactly as done for the LS case. The interaction parameters are 
again obtained from Goldschmidt. 35 The results are plotted in Figure 6. An attractive 
instability is found for J=4, a repulsive instability for J=2,0. The instability occurs at 
a value (Eq + E\ + 12E , 2)xo=l- This result is important in that it makes a very defi- 
nite prediction, if a paramagnon like picture analogous to He applies to heavy fermion 
superconductors, a pair state of relative J=4 should be realized. 

How does this change if we replace uranium by cerium? First, the energy zero does 
not contain Eq so that the bare interaction is more repulsive. Second, Eq is about 3 times 
larger since U is around 6 eV for cerium ions. 43 Although this means that the instability 
occurs for a smaller value of xo, we expect the interaction to be more repulsive since Eq 
is larger. This is illustrated in Figure 7 which is analogous to Figure 6 except parameters 
tabulated by Goldschmidt for the cerium ion are used 44 (with an Eq of 6.0 eV). As can 
be seen, the interaction in all channels is repulsive except very close to the instability for 
3—4. On the other hand, in strong coupling calculations for the paramagnon model in 
He, the calculated superfluid transition temperature actually turns off as the instability 
is approached since the energy scale of the paramagnon (proportional to 1 — Ixo) is going to 
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zero (that is, the T c maximum is near but not at the instability). 45 Because of this, pairing 
is possible for cerium alloys but much less likely than the uranium case where one finds a 
larger range of xo where there is attraction. An alternate view is that a pair wavefunction 
with two f quasiparticles has much less overlap with the bare f ion wavefunction in the 
cerium case since the / 2 component in cerium is much smaller than in uranium. This is 
in accord with experimental observations discussed in the first section. 

We now discuss the issue of orbital and momentum dependence of the susceptibility 
bubble. In heavy fermion uranium alloys, there is not much evidence for crystal field effects. 
This indicates that all the j=5/2 f orbitals are strongly mixed, as predicted by band theory 
(as discussed for Kondo lattice models by Zwicknagl, 46 if the Kondo temperature is larger 
than the crystal field splittings, then all f orbital energies get renormalized to the Fermi 
energy and are intermixed; this appears to be the case in U Pt%). Because of this, one would 
expect the orbital and momentum dependence of xo to be rather weak. There are inter- 
site interactions, though, which give the full susceptibility, x, momentum dependence. An 
argument against the importance of this effect for the pairing has been given by Anderson 2 
where he emphasizes the dominance of the on-site interaction given the large ratio of T c to 
Ep. The issue of intersite pairing will be discussed in the last section. On the other hand, 
the susceptibility in metals like UPt% is strongly dependent on field direction. Whether 
this is a consequence of crystal field effects or simply due to inter-site correlations is an 
unresolved matter 3 although neutron scattering data in UPts point to the latter. 6 If this 
is a crystal field effect, it can be represented by xo- If X is maximal for fields along the c 
axis (like in URv,2Si2), this indicates that the longitudinal bubbles are dominant. If one 
redoes Figure 6 with just longitudinal bubbles, then the J multiplets are split into various 
Mj terms (actually, those Mj combinations which have the appropriate crystal symmetry) 
with the maximum Mj configuration having the most attractive instability. On the other 
hand, if x is maximal for fields in the basal plane (like in UPt%) then transverse bubbles 
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with n=l are dominant (n as denned in Eq. 13 and 14). This in turn leads to the minimum 
Mj configuration being preferred. Higher values of n would indicate the importance of 
quadrupolar (n=2) effects, etc. These terms play an important role in certain theories. 14 ' 20 
These observations are summarized in Table 4 (where it should be noted that a J=4 state 
is always preferred). In this paper, these effects will be further ignored although they 
are relatively easy to incorporate if they can be accurately determined (and if they are 
indeed due to xo itself). An argument that these effects cannot be too strong is that 
an Mj=0 pair state is not consistent with experimental data in UPt% since it is a single 
dimensional group representation, although similar anisotropic spin fluctuation work in a 
spin-only approximation gave an analogous Mg=0 pair state, 10 which is consistent with 
observations of anisotropy in the upper critical field. 47 These issues will be discussed further 
in the next section. As will be discussed in the next section, projection of the vertex on 
the Fermi surface will also lead to lifting of the degeneracy of the J manifold (analogous 
to the lifting of L degeneracy in 3 ffe). It is that effect we concentrate on in this paper. 

We conclude this section with a discussion of the general vertex. In 3 ife, the full vertex 
can be written in a form analogous to the bare expression in Eq. 3, that is a density piece 
proportional to the delta functions and a spin piece proportional to the scalar product of 
spin operators. 27 Given the similarity of Eq. 3 to that of, say, Eq. 15, this should also 
be possible in the f electron case (as long as one restricts to states of definite LS or J). In 
particular, there will be a density piece, a quasi-spin (seniority) piece, and an orbital piece 
proportional to the scalar product of total angular momentum operators. This justifies 
some of the phenomeno logical interactions used in previous work. 13 Further exploitation of 
these ideas should give us a more fundamental insight into the properties of the full vertex 
for f electrons. We should note that the screened Slater integrals can be considered as 
analogues of the Landau F functions of Fermi liquid theory. 
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IV. Application to UPtj, 

We now wish to employ the formalism in the previous section to a real heavy fermion 
superconductor. We choose for this purpose UPt%. There are two good reasons for this. 
First, a variety of experimental data exist on this metal which gives us a fairly good idea 
about what the order parameter is. Second, extensive deHaas-vanAlphen data on U Pt%^ 
give a Fermi surface in fairly good agreement with LDA band structure calculations. 49 This 
indicates that the momentum dependence of the LDA wavefunctions is fairly trustworthy. 
The frequency dependence is not, of course, since the effective mass in the measurements 
is about 16 times the LDA band mass. This is due to the fact that in heavy fermions, 
the self energy has a large frequency derivative leading to a large mass enhancement. 
Investigation of transport properties indicates the momentum derivative of the self-energy 
must be rather weak so that the self-energy "rides" with the Fermi energy.^ That is why 
the shape of the LDA Fermi surface is about correct even though the mass is off by a 
large factor. These issues are of importance since we want to convert the formalism of 
the previous section to apply to quasiparticle states. We can do this approximately by 
taking the four bare external lines of the vertex in Eq. 1 and multiplying each of them 
by the wavefunction renormalization factor Z 1 / 2 , where as indicated above 1/Z~16. This 
represents the effect that only Z of the bare electron is in the quasiparticle pole. 

To proceed further, note that 

< J, a\r\f, a >= 5jj'5aa' r J,u (16) 

where a is a basis function of J which has the appropriate crystal symmetry. For axial 
symmetry, this would just be M. For hexagonal symmetry, they are listed in Table 5. 5 
The dependence of r on a occurs if anisotropy is put into xo as discussed in the last 
section. Also, for a multi-dimensional group representation, r will be the same for each 
a in the representation, unless the symmetry is lowered by some external perturbation. 
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This has relevance for Ginzburg-Landau models of the phase diagram for UPt% and will 
be discussed later. For now, though, we assume only a J dependence for r. These are 
plotted in Figure 6. Given this, we can now calculate the paring interaction on the Fermi 
surface. We do this by constructing the product \k,—k > and expanding this in terms 
of \J,a >. The first thing to note is that there are two degenerate states for each k (k, 
PTk where P is the parity operator and T the time reversal one) and two for —k (PA;, 
Tk). The combination ^(\k,Tk > — \PTk,Pk >) defines a pseudo-spin singlet, <1q. The 
three combinations \{\PTk,Tk > -\k, Pk >) = d x , -\{\PTk,Tk > +\k,Pk >) = d y , 
and \{\k, Tk > +\PTk, Pk >) = d z define a pseudo-spin triplet, known as the d vector. 

— * — * 

In the current theory, only the part of \k, —k > on the same site is involved in Eq. 16. 
(Note that although the pair interaction is only attractive for particles on the same site, 
the pairs are correlated out to a distance of the coherence length, much like the problem 
of bound states of a potential well where the particles spend most of their time outside 
the well. 29 ) Now, the part of \k > involving j=5/2 states is 52 

i^>=E<>>i ( 17 ) 

where \i runs from -5/2 to 5/2, i is the f atom site index (1,2 for UPt^), and n is the band 

index (five bands are predicted to cross the Fermi energy in UPt%). Thus, the coefficient 
— * — * 

of \k,—k > involving j=5/2 states on the same site with the correct group representation 
structure for a particular total J (J=0,2,4) is 

A t j = p ^Y.4<f k ( 18 ) 

where j represents the pseudo-spin combination (0 for singlet, x,y,z for triplet) and P is 
a projection operator which takes that part of the sum which has the form of one of the 
basis functions (Table 5) with the appropriate pseudo-spin combination discussed above. 
Because of antisymmetry, A changes sign from one site to the other site in the unit cell for 
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pseudo-spin triplets, and does not for pseudo-spin singlets (for one f atom per cell, only 
pseudo-spin singlets exist 53 ). Summarizing, the particle-particle vertex is 

< k', -k'\r\k, -k >= z 2 V rjA*./ aj, Ai aj (19) 

k 

Since Eq. 19 is separable in k and /c', this allows us to write down the BCS coupling 
constant 

\ Ja = NrjZ 2 Y,<\ A t ] \ 2 >k (20) 

3 

where N is the density of states, <>£ is an average over a narrow energy shell about the 
Fermi energy, and j runs over for even parity, x,y,z for odd parity. For U Pt$, this average 
was done on a regular grid of 561 k points in the irreducible wedge (1/24) of the Brillouin 
zone, keeping those nk states within 1 mRy of the Fermi energy (182 states total). Those 
points that are in symmetry planes of the zone are plotted in Figure 8. The number 1 
mRy was chosen so as to have enough points to give a good representation of the Fermi 
surface with this size grid. Note that in this model 

A Jaj (k)^A J ^ (21) 

where A is the order parameter. 

In Table 6, coupling constants 54 for UPt^ are shown modulo NTjZ 2 with the largest 
occuring for J=0, . This is just the BW state with spin-orbit, 55 with a coupling constant 
proportional to the square of the ratio of the j=5/2 f density of states to the total density 
of states. From Figure 6, though, the \2 a and Ao a coupling constants are repulsive and 
so do not play a role. For the attractive J=4 case, the largest coupling constants are for 
odd parity states. This is because three pseudo-spin triplet terms contribute to Eq. 20 
in this case as opposed to one pseudo-spin singlet for even parity states. The importance 
of this is that pseudo-spin triplets only exist because of the two f atoms per unit cell, 
which as illustrated in Table 1, all heavy fermion superconductors have. This is a property 



21 

of an on-site pairing theory, and has no relevance in near neighbor pairing models. The 
largest coupling constant occurs for i~jr (E2 U ) symmetry, although several other states 
have comparable sized coupling constants (r^,T^). Note that this is an odd parity two 
dimensional group representation. 

Let us discuss this state. For the odd parity case, the order parameter is a vector. 

— * 

At a general k, all three components are involved because of the relativistic nature of the 
wavefunction coefficients. For two dimensional group representations, this means that the 
state will be non- unitary (d x d* ^ 0). A discussion of the case of non- unitary d vectors 
can be found in Sigrist and Ueda 56 which play a major role in a recent phenomenological 
theory of UPt%. In the non-unitary case, there are two gaps for each k 

A a (k) 2 = \d(k)\ 2 + a\d(k) x d*(k)\ 2 (22) 

where a is +/-. Plots of this for the state are shown in Figures 9 and 10 with the 
total density of states shown in Figure 11. No attempt has been made to fit these gaps 
to simple functions 58 because the complex momentum dependence of the wavefunctions 
would prohibit this. Instead, some general properties can be inferred. For instance, the d z 
component vanishes for k z =0 as expected from earlier work. 59 A surprise, though, is that 
all three d vector components vanish on the zone face k z =iz/c. Thus, this state has a line 
node gap function, and provides a counterargument to earlier statements that line node 
gap functions are not possible for odd parity states. 59 Moreover, all d vector components 
vanish along the axis k x =0,k y =0 which gives rise to point nodes. A gap function of this 
sort (line nodes perpendicular to the c axis, point nodes along the c axis) is consistent with 
a variety of experimental data in UPt%, including specific heat, 60 transverse ultrasound, 61 
penetration depth, 62 thermal conductivity, 63 NMR, 64 and tunneling 65 data. Note that 
despite the small value of the second gap in Figure 10, no "normal" component is seen 
in the density of states of Figure 11. Thus, this non-unitary state differs from the one 
considered by Machida et al 57 where one of the gap components vanishes identically for 
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all k so that there is a normal component with half the value of the normal state (in the 
theory of Coleman et al, a similar normal component occurs). It should be remarked that 
although earlier specific heat data 60 indicated a sizable normal component, this is probably 
due to impurity effects since newer data do not show this component. 66 As for the two 
dimensional nature of the group representation, our current understanding of the field- 
temperature-pressure phase diagram of UPt^ is in strong support of such a state. 12 ' 67 ' 68 
In particular, weak magnetism is present which lowers the symmetry to orthorhombic. 
This acts to split the superconducting transition into two transitions. Pressure acts to 
eliminate both the magnetism and the splitting, 69 thus giving strong support for a two 
dimensional group representation, as opposed to two nearly degenerate single dimensional 
group representations. The main problem with this scenario is the presence of a term 
in the gradient part of the free energy which tends to mix the two components of the 
representation except for certain field directions. This is in contradiction to experiment 
which shows a fairly isotropic phase diagram with respect to field direction (this was the 
main motivation for the nearly degenerate model 70 ' 71 ). Sauls, though, has shown that 
for an axially symmetric Fermi surface and axially symmetric basis functions, this mixing 
term is zero for the E<i representation. 12 In the general case, it is not, but its mixing 
term is smaller than that of the E\ representation. This work will be discussed in another 
paper, 72 but it suffices to say here that for the current theory, the mixing term for our 
r§ (i?2u) state is small enough so as to be promising in regards to explaining the phase 
diagram. Alternate theories based on two nearly degenerate representations 70 ' 71 ' 73 are also 
consistent with the current theory given the closeness of the coupling constants for i~jT, 
i^ - , and r± (the latter state having the same line node structure as -Tg )• 

A final issue concerns the question of parity. No change below T c for UPt% has been 
seen in the Knight shift 74 or induced moment form factor, 7 ^ indicating no change in the 
magnetic susceptibility below T c . This is in support of an odd parity state, although 
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one could argue that in the heavy fermions, the quasiparticle (intraband) part of the 
susceptibility is small compared to the VanVleck (interband) part, so this conclusion is 
not definitive. Choi and Sauls 47 have also shown that the observed low temperature 
directional anisotropy of the upper critical field 76 is most easily explained with an odd 
parity pair state with Mg=0 (note, the quantization axis is assumed to be along c). Such 
a state came out of earlier non-relativistic spin fluctuation calculations which took into 
account the directional anisotropy of the susceptibility. 10 In the current case, though, spin 
and orbital components are mixed and so a pure Mg=0 state is not possible. On the 
other hand, the state found here, |Mj|=l, has the largest projection of J on the basal 
plane of any two dimensional group representation (Mj=0 has the largest projection, but 
is a singlet), so is promising in that regard. To test this quantitatively would require a 
detailed calculation of the upper critical field with both spin and orbital degrees of freedom 
taken into account. Certainly, if the current predicted state is correct, the H c 2 anisotropy 
cannot just be a spin effect as envisioned by Choi and Sauls. This can be seen as follows. 
A direct translation of their idea to the current state would be to have a state of pure d z 
character in pseudo-spin space. This can be tested by rotating in pseudo-spin space at 

— * — * 

each k so that the state \k > has maximal J z along the chosen quantization axis. This 
was done for quantization axes along the a, b, and c axes of the hexagonal crystal (for 
the a and b case, this lowers the system to orthorhombic symmetry). In all cases, the 
averages < |c^| 2 were within 20% to 30% of each other, i.e. there is no dominant 
d z component. This is consistent with the highly non-unitary nature of this state seen 
in Figures 9 and 10. As indicated above, then, orbital effects in the upper critical field 
must be playing an important role. This would require a theory for calculating H c 2 for a 
non-unitary multi-component d vector in the strong spin-orbit coupling limit. 

We finally turn to a discussion of T c . We should note that the density of states already 
contains the renormalization factor 1/Z (i.e., N = Nq/Z where Ng is the bare density 
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of states). This means that the prefactor NTjZ 2 in Eq. 20 is analogous to the form 
(1 + \z) found in strong coupling theories 45 since A^q-Tj would be the pairing coupling 
constant, A^, and 1/Z - 1 would be the mass coupling constant, Xz- This is consistent 
with the fact that renormalization of the external lines of the vertex is a strong coupling 
effect. On the other hand, the frequency dependence of the bubble has not been kept, so 
one has to simulate this by providing an energy cut-off of order the renormalized Fermi 
energy. We note that the size of the specific heat coefficient 3 and the neutron scattering 
linewidth^ are consistent with a renormalized energy scale for UPt% of order 5 meV. This 
is also consistent with estimating a band structure Fermi energy and multiplying this by Z. 
The last thing to be determined is This is made difficult by the fact that xo is being 
treated as a number in this paper, whereas in reality it is a higly complicated function 
of momentum, frequency, and band and orbital indices. Given the strong dependence 
of Fj on xo (and also the question of the energy zero), the exponential dependence of 
T c on coupling constant, and the uncertainties mentioned in the above approximations, 
the most illustrative approach is just to see what value of xo is needed to obtain the 
observed T c . With a T c of 0.5 K and a cut-off of 5 meV, this gives a value for A of 0.205. 
Since A = NT A Z 2 c where c from Table 6 for J=4, r Q is 0.125, this gives a T 4 of -2.1 
eV relative to the energy zero. From Figure 6, this gives a xo °f 0.335, which is 0.92 
of the divergence value, giving a Stoner renormalization of 12, comparable to the mass 
renormalization value of 16 assumed from the beginning. Since spin fluctuation models 
based on the observed heavy fermion dynamic susceptibility have a mass renormalization 
which goes like the Stoner factor, as opposed to the log of the Stoner factor which one 
gets for a Lindhard function, 8 then there is a overall consistency in these numbers. This 
can be further demonstrated by estimating x, obtained by multiplying xo by the Stoner 
renormalization (12) then by the square of the orbital degeneracy (36). This gives a value 
of about 150 states/eV, comparable to the 180 states/eV given by the specific heat 7 of 
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420 mJ/mol K 1 ? Moreover, converting x to proper units (with g—6/7 and j=5/2) gives 
0.0075 emu/mol, comparable to data from susceptibility measurements. 3 To obtain a more 
reliable estimate of these parameters would require doing a full strong coupling calculation 
retaining the frequency dependence of the bubble. 45 

The large estimated size of Ij of order 2 eV which must be renormalized downwards by 
Z 1 might seem somewhat worrisome. After all, wouldn't one expect high T c in transition 
metals where Z is closer to one? The question of renormalizing the interaction downwards 
has been discussed by Anderson 2 and reviewed by Lee et al. 77 The main point to emphasize 
here is that the interaction parameters of this paper are only appropriate for a system close 
to an / 2 configuration, and the same thing that is causing the large value of i~j is also 
causing a small value for Z. In transition metals, the Ej~ parameters are largely screened 
out and play no role. Instead, one collapses back to a standard spin fluctuation model with 
a Stoner interaction parameter I. Estimates based on this I give extremely low estimates 
of T c , even in palladium which has a large Stoner renormalization. 78 

V. Future Directions 

An advantage of the current approach is that the theory is systematically improvable 
by removing various approximations made in this paper. The most severe of these is 
treating xo as a number. A proper strong coupling calculation would include the frequency 
dependence of this object. This is not too difficult if the simple relaxational form is used 

XOH = (23) 

where F is the neutron scattering linewidth.^ Work of this sort has been done in earlier spin 
fluctuation models. 8 ' 10 Of more interest is the momentum and orbital dependence of this 
object. The philosophy of this paper is similar to that espoused early on by Anderson, 2 that 
is the size of the effective interaction is large enough (large ratio of T c to Ep) that the on- 
site interaction must play the central role. This is also consistent with the observation that 
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one of the denning properties of the heavy fermion metals is their large f atom separation 79 
giving weak inter-site interaction effects. Even neutron scattering data indicate that that 
part of the susceptibility where the bulk of the fluctuating moment is has a relatively 
mild momentum dependence. 6 On the other hand, these effects played a crucial role in 
earlier spin fluctuation models, so it is of interest to see how they would enter the current 
formalism. In the real space approach taken here, these effects would be simulated by 
adding a term xo(R, R') m addition to the term xo(Rj R) where R, R' are site indices with 
R 1 a near neighbor of R. Near neighbor effects would show up in Eq. 1 at first order in 
this bubble. At second order in this bubble, there would be terms which would also affect 
on-site pairing (xo(R, R')xo(R'i R))- Assuming Eq. 1 can be solved, one is left with a 
vertex which contains on-site terms, near neighbor terms, etc. The general properties of 
this vertex along with the functional forms of the near neighbor and next near neighbor 
pairs are discussed in Appel and Hertel. 5 One simply has to combine this work with that 
one to get a complete solution. It is complicated, but doable. A simple argument, though, 
can be used against near neighbor pairing in that the same sign difference between on- 
site and near neighbor interactions in the antiferromagnetic case which gives rise to S=0 
pairing in the single orbital case (end of Section II) should give J=0 pairing in the orbitally 
degenerate case, since from Figure 6, J=0 has maximal repulsion for the on-site case. Such 
a state is a singlet, and is completely inconsistent with the data we have discussed on 
UPt 3 . 

The other issue concerns orbital dependence of the bubble. Calculations of this sort 
exist in the literature for U Pt^° and are rather tedious, as they involve calculating matrix 
elements of relativistic wavefunctions over a fine enough grid in k space to get reliable values 
of xo- As discussed in Section III, one might get around this difficulty by simulating these 
effects with some effective crystal field model. The most probable picture based on the 
temperature dependence of the neutron scattering data, though, is that this effect enters 
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most prominently in the xo(R,R') term. 6 

The next issue concerns feedback effects which are important in the physics of super- 
fluid ^He. In that case, all L=l states are degenerate at T c . Therefore, one would expect 
the isotropic state, the BW state, to have the lowest free energy. On the other hand, the 
susceptibility changes below T c which affects the pair interaction and favors states with 
maximal anisotropy, giving rise to the ABM state. 27 ' 29 In the current problem, though, (1) 
this degeneracy is already broken in the normal state due to the crystal lattice and (2) there 
is no experimental evidence for a change in the susceptibility below T c (based on Knight 
shift 74 and induced moment form factor 75 measurements), although it should be noted 
that the interband, or VanVleck, component most likely dominates the susceptibility. 80 
(Neutron scattering experiments which access the low momentum, low frequency part of 
the dynamic susceptibility would be helpful in extracting out the quasiparticle part of 
the susceptibility and seeing how it changes below T c .) Because of this, feedback effects 
probably do not play an important role in the heavy fermion problem (the large ratio of 
T c to Ep would argue against this, though). A related effect is whether the interaction 
changes as a function of field (this could be connected to the H c 2 anisotropy discussed in 
the previous section). Magnetization data look very linear in field for all field directions 
for the field range of interest 3 which would argue against this. On the other hand, as T c 
depends exponentially on coupling constant, small changes in the quasiparticle wavefunc- 
tions in an applied field could lead to noticable effects, especially given the low effective 
Fermi energy. This could be simulated in the current theory by rediagonalzing the band 
structure wavefunctions in the presence of the appropriate field and see how the coupling 
constants listed in Table 6 change. 

A related effect is the observed splitting of T c in UPt% that was discussed in the 
previous section. This has been treated in the past by a phenomenological symmetry 
breaking field thought to be due to the orthorhombic distortion associated with the weak 
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antiferromagnetism. 81 This could also be simulated in the current model by applying a 
weak staggered field (the new orthorhombic cell would contain four uranium atoms) and 
rediagonalizing the band structure wavefunctions, from which the splitting of the coupling 
constants could be determined. 

The author would like to conclude by saying that the physics of heavy fermion super- 
conductors is complicated enough (as this paper demonstrates) that the picture offered 
here may not be complete. On the other hand, he feels that the ultimate theory for these 
materials must look in some form like what is being proposed here, since the orbital de- 
pendence of the f electrons and the short range nature of the interactions should play a 
crucial role. It is promosing that this model has certain qualitative features reflected in the 
data (preference for uranium with two f atoms per unit cell) which are hard to understand 
from earlier spin fluctuation theories. Also, the predicted pair state for UPt% has many 
promising aspects also missing in earlier theories. Moreover, there is a conceptual beauty 
to having a pair state which has maximal L and maximal S, as this is a direct general- 
ization of the physics of ^He with which heavy fermions share many qualitative features. 
Hopefully, with increased experimental and theoretical effort, we can determine whether 
this is indeed the right approach to pursue for solving this problem. 
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Table 1. List of known heavy fermion superconductors with the number of f atoms per 
unit cell. In parenthesis is the nature of the low temperature distorted phase in the single 
f atom case (QP - quadrupolar, AF - antiferromagnetic, ? - not fully determined) and the 
resulting number of f atoms. 



Case 


f Atoms 


UPt s 


2 


UBe 13 


2 


U 2 PtC 2 


2 


URu 2 Si 2 


1 (QP/AF - 2) 


UPd 2 Ah 


1 (AF - 2) 


UNi 2 Al 3 


1 (AF - 2) 


CeCu 2 Si 2 


1 (? - 2 ?) 
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Table 2. / 2 energies in the LS scheme. 37 



Term 


Energy 




Eq - 9E 3 


3 F 


Eq 


3p 


Eq + 33£ 3 


l I 


E + 2Ei + 70E 2 + 7E 3 


*G 


E + 2Ei - 260E 2 - 4E 3 




E + 2Ei + 286E 2 - 11E 3 




E + 9Ei 



37 

Table 3. / 2 energies in the JJ scheme. 42 Note that the parameters are different from 
those defined in the LS scheme. 

Term Energy 
J = 4 E -5E 2 
J = 2 E + 9E 2 

J = ^0 + 3^i 
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Table 4. Summary of anisotropic r for J=4 using parameters of Figure 6. n is the type of 
bubble used (defined in Eq. 14), xo is the value (eV~ l ) at which the divergence occurs, 
and M signifies which Mj state is the most attractive (with other attractive Mj states 
listed in parenthesis). Note there is no attraction for n=5. For J=2, attraction is found 
for Mj = in the n=3,4 cases, but weaker than J=4. For J=0, no attraction is found. 



n 


XO 


M 





0.365 


4(3) 


1 


0.365 


(1,2) 


2 


0.420 


0(2) 


3 


0.420 


1 


4 


0.420 





5 


0.420 
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Table 5. Hexagonal basis functions for J=4. The forms listed in this table should be 
(a) antisymmetrized > \v > —\v > |/z >) and (b) symmetrized (+ representation) 
or antisymmetrized (- representation) with respect to site before use. For A, a and f3 
are variational coefficients such that the sum of their squares is equal to one, and this 
representation occurs twice (a, (3 and (3, —a). Note that A and A are doublets obtained 



by replacing \/i > by | — /i >. 

Rep Basis Function 

r 5 a\5/2 > |3/2 > +/3(0.8018|5/2 > | - 1/2 > +0.5976|3/2 > 1 1/2 >) 

r 3 0.7071|5/2 > 1 1/2 > +0.7071| - 5/2 > | - 1/2 > 

A 0.7071|5/2 > |l/2 > -0.7071| - 5/2 > | - 1/2 > 

r 6 0.5345|5/2 > I - 3/2 > +0.8452|3/2 > | - 1/2 > 

A 0.2673(|5/2 > I - 5/2 > +3|3/2 > | - 3/2 > +2|l/2 > | - 1/2 >) 
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Table 6. Coupling constants for UPt%. These should be multiplied by the quantity NT^Z 
to convert to real coupling constants. 



Rep 


J=4 (+) 


J=4 (-) 


J=2 (+) 


J=2 (-) 


J=0 (+) 


J=0 (-) 


A 


0.069 


0.073 










A 


0.029 


0.101 


0.049 


0.071 






A 


0.024 


0.064 










A 


0.013 


0.120 










A 


0.018 


0.125 


0.035 


0.106 






A 


0.076 


0.114 


0.065 


0.099 


0.495 


0.057 
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FIGURE CAPTIONS 



1. Interaction potential between two He atoms. The interaction potential of two f 
electrons of a uranium ion would look similar with appropriately scaled axes (with 
attraction in that case due to the ion core). 

2. f 2 multiplet structure of a £/" 4+ ion. 35 ' 82 Energy is plotted versus J, with labels 
referring to L. 

3. Effective interaction of Eq. 10 versus xo- 

4. Effective interaction of Eq. 11 (triplet, lower curve) and Eq. 12 (singlet, upper curve) 
versus xo- This is very similar to the effective interaction in the 3 iiZe problem. 

5. Effective interaction (LS) in eV for 3 i7, 3 F, 3 P, and 1 7 versus xo f° r parameters 
appropriate to a U ion 35 ' 40 (E Q = 1225 meV, E 1 = 470.3 meV, E 2 = 1.923 meV, E 3 
= 43.28 meV). The zeros of energy for the f 1 and f 2 cases are marked by the dashed 
lines. 

6. Effective interaction (JJ) in eV for J=4,2,0 versus xo f° r a U ion. Same parameters 
and notation as in Figure 5. 

7. Effective interaction (JJ) in eV for J=4,2,0 versus xo f° r parameters appropriate to 
a Ce ion 44 ' 43 (Eq = 6000 meV, E\ = 484.5 meV, E 2 = 2.293 meV, E 3 = 47.67 meV; 
note these are LS E^). The zero of energy is the dashed line. 

8. Plot of the k points used in the calculations on UPt% in symmetry planes of the 
Brillouin zone constructed from a uniform grid within an energy cutoff of 1 mRy of 
the Fermi energy. Lines represent the Fermi surface. 
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9. Plots of |Aj(fc)| with a = + for r Q state on the grid of points of Figure 8. Size of 
the dots represent the magnitude of the gap. Where no dots appear, gap is zero or 
very small. 

10. Same as Figure 9, but with a = -. 

11. Smoothed fit to the density of states (normalized to the normal state value) con- 
structed from the gaps plotted in Figures 9 and 10. Energy units are normalized to 
the value of the largest gap. 
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